SMOOTH ATTRACTORS FOR THE BRINKMAN-FORCHHEIMER 
EQUATIONS WITH FAST GROWING NONLINEARITIES 
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Abstract. We prove the existence of regular dissipative solutions and global attrac- 
tors for the 3D Brinkmann-Forchheimer equations with the nonlinearity of an arbitrary 
polynomial growth rate. In order to obtain this result, we prove the maximal regularity 
estimate for the corresponding semi-linear stationary Stokes problem using some mod- 
ification of the nonlinear localization technique. The applications of our results to the 
Brinkmann-Forchheimer equation with the Navier-Stokes inertial term are also consid- 
ered. 
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We study the Brinkman-Forchheimer (BF) equations in the following form: 



1. Introduction 



d t u — Am + f(u) + Vp = g, divu = 0, 
u\ dn = 0, u\ t=Q =u . 

Here Q C M 3 is an open, bounded domain with C 2 boundary dQ, g = g(x) = (<7i, <?2, #3) 
is a given function, u = w 2 , is the fluid velocity vector, p is the pressure and / is 
a given nonlinearity. 

The BF equations are used to describe the fluid flow in a saturated porous media, see 
[T6l [21] and references therein. The typical example for / is the following one: 

f(u) = au + 6|u| r_1 u, r G [1, oo), (1.2) 

where a G R and b > are the Darcy and Forchheimer coefficients respectively (the 
original Brinkman-Forchheimer model corresponds to the choice r = 2, more complicated 
nonlinear terms (r ^ 2) appear, e.g., in the theory of non-Newtonian fluids, see [20]). Note 
also that the analogous equations are used in the study of tidal dynamics (see [S].[T3]). 

Number of papers is devoted to the mathematical study of of the BF equations, for 
instance, continuous dependence on changes in Brinkman and Forchheimer coefficients 
and convergence of solutions of BF equations to the solution of the Forchheimer equation 

d t u + f(u) + Vp = g, div u — 0, 

as the viscosity tends to zero have been established in [2J, [31 [12l HU [21] (see also references 
therein), and the long-time behavior of solutions for (11.11) has been studied in terms of 
global attractors in [17] [23] and [23]. However, to the best of our knowledge, only the 
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case of the so-called subcritical growth rate of the nonlinearity / (r < 3 in (11.2jl ) has been 
considered in the literature. 

The main aim of the present paper is to remove this growth restriction and verify the 
global existence, uniqueness and dissipativity of smooth solutions of the BF equations for 
the large class of nonlinearites / of the arbitrary growth exponent r > 1. 

Namely, we assume that / G C 2 (IR 3 ,IR 3 ) satisfies the following conditions: 

1) f (u)v.v > (-K + K\u\ r ~ 1 )\v\ 2 , \/u,v G R 3 , 

2) \f{u)\ < C(l + l^r" 1 ), VmgM 3 , ^ ' ' 

where K, C, k are some positive constants, r > 1 and u.v stands for the standard inner 
product in M 3 . 

Our key technical tool is the maximal regularity result for the stationary problem 

- Aw + f(w) + Vp = g, divw = 0, u\ m = (1.4) 

which claims that the solution w belongs to H 2 if g G L 2 . This result is straightforward 
for the case of periodic boundary conditions (it follows via the multiplication of the 
equation by Aw and integrating by parts). However, for the case of Dirichlet boundary 
conditions it is far from being immediate since the additional uncontrollable boundary 
terms arise after the multiplication of the equation by Aw and integrating by parts. 
Following the approach developed in [9], we overcome this problem using some kind of 
nonlinear localization technique, see Apendix below. 

In addition, we apply our maximal regularity result in order to establish the existence 
of smooth solutions for the so-called convective BF equations: 

d t u + (u, V)u — Au + f(u) + Vp = g, divu = 0, . _ N 

i i (1-5) 

u \dn = °> M lt=o = u ° 

under the assumption f ll.2p with r > 3. Note that the case / = corresponds to the clas- 
sical Navier-Stokes problem where the existence of smooth solutions is an open problem. 
However, as also known (see [HJ]) the sufficiently strong nonlinearity / produces some 
kind of regularizing effect. Again, in contrast to the previous works, no upper bounds for 
the exponent r are posed here. 

The paper is organized as follows. A number of a priori estimates which are necessary 
to handle equation (11 .ip is given in Section [2j Existence, uniqueness and regularity of 
solutions for the BF equation as well as the existence of the associated global attractor 
are established in Section These results are extended to the case of convective BF equa- 
tions (jl.5p in Section |H Finally, the crucial maximal regularity results for the stationary 
equations (II. ip and (II. 5p are obtained in Appendix. 

2. A PRIORI ESTIMATES. 

In this section, we obtain a number of a priori estimates for the solutions of the prob- 
lem (II. ip assuming that the sufficiently regular solution (u,p) of this equation is given. 
These estimates will be used in the next sections in order to establish the existence and 
uniqueness of solution, their regularity, etc. 
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We start with introducing the standard notations. As usual, we denote by W l ' p (Q) 
the Sobolev space of all functions whose distributional derivatives up to order I belong to 
LP (ft). The Hilbert spaces W l > 2 {tt) will be also denoted by Lf'(O). 

For the vector valued functions v = (v i, t>2, vs), and u = (ui,U2, 1*3) we denote by (u, v) 
the standard inner product in [L 2 (ft)] 3 : 

3 

(u,v) := J2( v i> u i"> L2 ( u )> 

3=1 

3 

and write ||Vw|| 2 2 instead of ll^^llia- In the sequel, where it does not lead to mis- 

understandings, we will also use the notation H l (Q) and W l,p (Q) for the spaces of vector 
valued functions [if' (ft)] 3 and [W z ' p (ft)] 3 respectively. 
As usual, we set 

V := {v G (C™(n)) 3 : divv = 0} , 

and denote by H and H\ — V the closure of V in L2(ft) and if 1 (ft) topology respectively. 
And, more generally, H s := D(A S ^ 2 ), where A := IIA and II is the classical Helmholz- 
Leray orthogonal projection in L 2 (ft) onto the space H. In particular, since Q is smooth 
and bounded, we have 

H = {u G L 2 (ft), divu = 0, (u,n)\ dn = 0}, H 1 := H n H%(Q), H 2 = H 1 n H 2 (£l), 

see e.g. [10J. 

The next lemma gives the usual energy estimate for the BF equation. 

Lemma 2.1. Let (u,p) be a sufficiently smooth solution of problem (11.11) . Then the 
following estimate holds: 

\\u(t)\\ 2 L2 + J^ +1 [\\Vu{s)f L , + \\u(s)\\ r L tl] ds < C\\u(0)\\he- at + C(l + \\g\\h), (2.1) 

where the positive constants C and a are independent of t and the concrete choice of the 
solution (u,p). 

Proof. Indeed, multiplying equation (11.11) by u, integrating over x G ft, using that 
f{u).u > —C + k|-u| t ' +1 and (Vp, u) = (p, divit) = and arguing in a standard way, 
we have 

\d t \\u{t)\\ 2 L2 + a\\u(t)\\m + "KOKi < C{1 + \\g\\l>) (2.2) 

for some positive a and C which are independent of u and t. Applying the Gronwall 
inequality to the last estimate, we derive (12.11) and finish the proof of the lemma. □ 

Remark 2.2. The standard (for the reaction-diffusion equations) next step in a priori 
estimates would be the multiplication of equation (II .ip by Au (or tAu) and obtaining 
the dissipative estimate in H l together with the L 2 — > H l parabolic smoothing property. 
However, in our case, this scheme looks not applicable since Att| an 7^ in general and the 
term with pressure will not disappear. Multiplication by TIAu (where II is the Helmholz- 
Leray projector to the divergent free vector fields) also does not work due to the presence of 
the non-linearity / with arbitrary growth rate. So, we have to skip this step and estimate 
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the L 2 -norm of d t u instead differentiating equation by t and using the quasi-monotonicity 
of /. The H l (and H 2 ) estimate will be obtained after that using the maximal regularity 
theorem for the elliptic problem ( 15. II) . see Appendix). 

The next simple corollary is, however, crucial for our method of proving the existence 
and dissipativity of the iY 2 -solutions. 

Corollary 2.3. Let (u,p) be a sufficiently regular solution of problem ( II. lj) . Then, the 
following estimate holds: 

\\dtu\\ L , m+llH -*) < Q(\\u(0)\\ L 2)e- at + Q(\\g\\ L2 ), (2.3) 

where the monotone function Q and the constant C are independent of t and u. 

Proof. Indeed, applying the Helmholz-Leray projector II to both sides of equation (11.11) 
and using that div dtu = 0, we arrive at 

d t u = Au- Uf(u) + Tig. (2.4) 

Thanks to the growth restriction on / and the control (12. ip . we have 

ll/WHWw-i],^) < C\\u(0)\\l 2 e- at + C(l + \\g\\l 2 ) 

with r* := ^-y-. Using now that the Helmholz-Leray projector II : U* — > U* together 
with the embedding U C H~ 2 (recall that n = 3), we arrive at 

l|n/(u)|Ui ([titfl]iH - 3) <Q(\H0)\\L*)e~ at + Q(\\g\\Li) 

for some monotone increasing function Q. This estimate, together with (12.41) and the 
control of u given by the energy estimate (12.11) give the desired estimate (I2.3P and finish 
the proof of the corollary. □ 

Let us now differentiate (11.11) with respect to time and denote v = d t u. Then, this 
function solves 

d t v = Av -f(u)v + Vq, divv = 0, v(0) = Au(0) -Uf(u(0)) + Ug. (2.5) 
Moreover, using the embedding H 2 C C, we see that 

\\v(0)\\ L * <Q(\H0)\\H*) + \\g\\v (2.6) 

and, therefore, the L 2 -norm of the initial data for v is under the control if w(0) e H 2 . 
The next Lemma gives the control of v(t) for all t > 0. 

Lemma 2.4. Let (u,p) be a sufficiently regular solution of problem (11.11) . Then, the 
following estimate holds: 

rt+l 

\\v(t)\\h + J \Hs)\\^ds<Q(\\um H >)e Kt + Q(\\9\\l^ (2.7) 
for some positive constant K and monotone function Q. 

Proof. Multiplying equation ( 12. 5 p by v(t), integrating over fl and using that (f'(u)v) -v > 
— K\v | 2 , Wu, v G M 3 (see the condition (11.31) ). we arrive at 

dt\\vmh + Mt)\\m<2K\\v(t)\\ 2 L2 . (2.8) 
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Applying the Gronwall inequality to this estimate, we arrive at (I2.7P and finish the proof 
of the lemma. □ 

Corollary 2.5. Let (u,p) be a sufficiently smooth solution of the problem (11. ip . Then, 
the following estimate holds: 

\\u{t)\\ H * + \\V P (t)\\ L , < Q(\\u(0)\\ H2 )e Kt + Q{\\g\\v) (2.9) 

for some positive constant K and monotone function Q independent oft and uq. 

Indeed, due to the control (12.71) . we may rewrite equation (11.11) as an elliptic boundary 
value problem 

Aw(t) - f(w(t)) + Vp(t) = g u (t) := -g + d t u(t) (2.10) 

and apply the maximal regularity result of Theorem 15.21 (see Appendix) to that equation. 
Together with (12.71) this gives indeed estimate (12.91) and proves the corollary. 

We, however, note that the proved estimate (I2.9P is divergent as t — > oo and, by that 
reason, is not sufficient to verify the dissipativity of the problem (11.11) in H 2 . In order to 
overcome this drawback, we need the L 2 — > H 2 smoothing property for the solutions of 
(jl.ip . This result will be obtained exploiting the parabolic smoothing for equation (12.51) 
together with the already established control (I2.3P for v (t) = dtu(t). 

Lemma 2.6. Let (u,p) be a sufficiently regular solution of the problem (11.11) . Then, the 
following estimate holds: 

\\d t u(t)\\ L 2 < i±£ (Q(|| M (0)|| L2 )e- Qi + Q(\\g\\v)) , t > 0, (2.11) 

where the positive constant a and the monotone function Q are independent of t and u. 

Proof. We first note that, due to the energy estimate (12. ip . it is sufficient to verify (12.111) 
for t e (0, 1] only. To this end, we multiply (12. 8ft by t (where the exponent will be 
specified later) and integrate with respect to t. Then, we have 

sup {s N \\v(s)\\ 2 L2 }+ [ t s N \\v(s)\\ 2 H1 ds<C /V-X^li, :=/(t), (2.12) 
se[o,t] Jo Jo 

where C = C(N, K) is independent of t and u. 

We estimate I(t) using (12.31) and the interpolation inequality ||f ||l 2 < CII^IIh 3 """" 



\n IP- 



lit) <C sup {s N/2 \\v{s)\\ L 2} [ s N/2 - x \\v{s)\\ L 2 ds< 
se[o,t] Jo 



se[o,t] 

< sup {s N / 2 \\v(s)\\ L i} A^ii^,)!!^) 273 ^ 72 - 3 ^^)!!^) 173 ^^ 

se[o,t] Jo 

< 1/2 sup {^11^)1114 + 1/2 f s N \\v(s)\\ 2 Hl ds+ 
se[o,t] Jo 

+ C ( t s N / 2 - 3 \\v(s)\\ H - 2 ds) . (2.13) 
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Fixing now N = Q, using the control ( 12 .3p in order to estimate the right-hand side of 
(I2.13P and inserting it into the right-hand side of (I2.12p . we see that 

sup {s 6 \\v(s)\\ 2 L ,} < 1/2 sup {s 6 \\v(s)\\ 2 L2 } + Q(\\u(0)\\ L 2) + Q(\\g\\ L2 ) (2.14) 

se[o,t] s€[0,t] 

(recall, we have assumed that t < 1). It only remains to note that (12.141) immediately 
gives (12.111) for t < 1. Lemma [2.61 is proved. □ 

We summarize the obtained estimates in the following theorem. 

Theorem 2.7. Let (u,p) be a sufficiently regular solution of the problem (11.11) . Then, 
the following estimate holds: 

\\u{t)\\ m + ||Vp(*)||jffi < Q(\\u(0)\\ m )e- at + Q(\\g\\ L 2), (2.15) 

where the positive constant a and a monotone function Q are independent oft and u. 
Moreover, the following smoothing property is valid: 

\\u(t)\\m + ||Vp(t)|| L2 < Q (i±^|| u (0)|U 2 ) e- at + Q(\\g\\ L2 ), t > 0. (2.16) 

Indeed, the estimate (I2.16P is an immediate corollary of ( 12. lip and the maximal elliptic 
regularity of Theorem 15.21 applied to the elliptic equation (I2.10p . In order to verify (I2.15p . 
it is sufficient to use the divergent in time estimate ( 12 .9p for t < 1 and estimate (I2.16P for 
t > 1. 

3. Well-posedness and attractors 

The estimates obtained in the previous section, allow us to prove the existence and 
uniqueness of a solution of the problem ( II. ip as well as to establish existence of the global 
attractor for the associated semigroup. We start with the definition of a weak solution of 
that equation excluding the pressure in a standard way. 

Definition 3.1. A function 

u e C([0, oo), H) n LlM oo), H 1 ) n Lff([0, oo), (3.1) 

is called a weak solution of (II. ip if it satisfies ( 12. 4 p in the sense of distributions, i.e., 

- / (u(t), d t cp(t)) dt = - I (Vu{t), Vcp(t)) - (/(«(*)), <p(t)) + {g, <p(t)) dt 
Jr Jr 

for all Lp G C£°(IR + x Q) such that div (p(t) = 0. 

The next lemma establishes the uniqueness of a weak solution. 

Lemma 3.2. Let the nonlinearity f satisfy assumptions (II .3p . Then, the weak solution 
of problem (II. ip is unique. Moreover, for any two solutions Ui(t) andu2(t) (with different 
initial data) of the equation (11.11) . the following estimate holds: 

H(f) -u 2 (t)\\ L 2 < e^- Al )*|M0) -n 2 (0)|| x2 , (3.2) 
where K is the same as in (11.31) and \\ > is the first eigenvalue of the operator A. 
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Proof. Let U\{t) and u 2 {t) be two different energy solutions of (II. II) and let v(t) := Ui(t) — 
U2(t). Then, this function solves: 

d t v = Av-U(f(u 1 )-fM) 1 v(0) = Ul (0)-u 2 (0). (3.3) 

Note that, due to the regularity (13. ip of a weak solution and the growth restrictions on 
/, all terms in equation (13.31) belong to the space 

L 2 ([0, T], H~ l ) + L 1+1/r ([0, T], L 1+1,r \9)) = [L 2 ([0, T], H 1 ) n L r+1 ([0, T], L 1+r (ft))]*. 

In particular, the function t — > \\u(t)\\ 2 H is absolutely continuous and 

j t \\u(t)\\l 2 = 2(d t u(t)Mt)). 
Multiplying now equation ( 13. 3ft by v(t), integrating over Q and using the inequality 

(f(ui) - f(u 2 (t)).( Ul - u 2 ) > -K\ Ul - u 2 \ 2 , Vwi, u 2 G M 3 
(due to the first assumption of ( 11.31) ). we arrive at 

V*j t \Ht)\\h < K\W)\\h ~ (Mt)Xt)) < IK- Ax)|k(t)||| 2 (3.4) 

and the Gronwall inequality now gives the uniqueness and estimate (I3.2p . Lemma [3.21 is 
proved. □ 

We are now able to state our main result on the well-posedness and regularity of solu- 
tions of problem (11.11) . 

Theorem 3.3. Let the nonlinearity f satisfy assumptions (11.31) and let g G L 2 (Q). Then, 
for every Uq G H, problem (II. ip possesses a unique weak solution u (in the sense of 
Definition (13. ip ). Moreover, u(t) G H 2 for all t > and the estimate (12.161) holds. In 
addition, if u G H 2 , the estimate (12.151) also holds. 

Proof. Indeed, the existence of a weak solution can be obtained in a standard way using, 
say, the Galerkin approximation method. The uniqueness is proved in Lemma [3721 Thus, 
we only need to justify the estimates (I2.16P and (I2.15p . To this end, we note that the 
estimates (12.71) and (12. lip for the differentiated equation (I2.5P can be also first obtained 
on the level of the Galerkin approximations and then justified by passing to the limit 
(remind that the uniqueness of a weak solution holds). Finally, rewriting the problem 
(II. ip in the form of elliptic problem (I2.10p and using the Theorem 15. 2| we justify the 
desired estimates (12.161) and (I2.15p . Thus, Theorem 13.31 is proved. □ 

Thus, under the assumptions of Theorem l3.3[ the Brinkman-Forchheimer problem (11.11) 
generates a dissipative semigroup S(t) in the phase space H: 

S(t) :H^H, S(t)u := u(t), (3.5) 

where u(t) solves (11.11) with u(0) = uq. Our next task is to verify the existence of a global 
attractor for that semigroup. For the convenience of the reader, we start with reminding 
the definition of the attractor, see [1] , [7] , [TTJ , [22] for more details. 
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Definition 3.4. A set A C H is a global attractor of a semigroup S(t) : H — y H if the 
following properties are satisfied: 

1) A is a compact subset of if; 

2) ^4 is strictly invariant: S(t)A = A for all t > 0; 

3) It attracts the images of all bounded sets as time goes to infinity, i.e., for every 
bounded subset B C H and every neighborhood 0(A) of A, there exists T = T(_B, 0) 
such that 

S(t)B c 0(„4), Vt > T. 

The following theorem states the existence of the attractor for the problem considered. 

Theorem 3.5. Let the assumptions of Theorem \3. 3\ hold. Then the solution semigroup 
( 13. 5 p associated with the Brinkman-Forchheimer equation ( II. ip possesses a global attractor 
A (in the sense of the above definition) which is bounded in H 2 and is generated by all 
complete bounded solutions of ( II. ip defined for all t e R; 

A = K\ t=0 , (3.6) 

w/iere /C := {u E C b (R,H 2 ), u solves (TTTTj) }. 

Indeed, according to the abstract attractor existence theorem (see e.g., [TJ,[22]), we 
only need to check that the considered semigroup is continuous with respect to the initial 
data (for every fixed t) and it possesses a compact absorbing set in H. But the first 
assertion is an immediate corollary of Lemma 13.21 and the second one follows from the 
estimate (I2.16p . Moreover, this estimate gives the absorbing set bounded in H 2 . Since 
the attractor is always contained in an absorbing set, we have verified the existence of 
a global attractor A which is bounded in H 2 . Finally, the representation ( 13. 6 j) of the 
attractor in terms of completer bounded trajectories is also a standard corollary of the 
attractor existence theorem mentioned above. 

Remark 3.6. Although, we have stated only the if 2 -regularity of the attractor A, it can 
be further improved (if f, Q and g are smooth enough) using the maximal regularity for 
the linear Stokes equation and bootstrapping. In particular, if /, Q and g are C°° smooth, 
the attractor will be also C°°-smooth. 

Another standard corollary of the general theory is the fact that the obtained attractor 
has a finite Hausdorff and fractal dimension in H. The proof of this fact is a straight- 
forward implementation of the volume contraction technique to our equation (see e.g., 
[TJ |22]). Indeed, due to the embedding H 2 C C, the nonlinearity / is subordinated to 
the linear part of the equation (no matter how large is the growth exponent r) and one 
even is able to reduce formally the problem considered to the case of abstract semilinear 
parabolic equations. 

To conclude this section, we discuss the particular case of (II. ip where 

f(u) = -V u F(u), (3.7) 

for some scalar function F G C 2 (M 3 ). Note that this condition is satisfied for the "most 
natural" nonlinearities 

f(u) = au\u\ — bu. 
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In that case, multiplying the equation by d t u and integrating over Q, we get 

j t £(u(t)) = -\\d t u(t)\\h < 0, 

where 

C(u) :=i(Vw,V«) + (F(«),l). 

Thus, the solution semigroup S(t) possesses the global Lyapunov functional C(u) and 
applying the standard arguments (see [7].|11]) to our problem, we obtain the following 
result. 

Corollary 3.7. Let the assumptions of Theorem \3.5\ and the condition (13. 7p be satisfied. 
Then, every trajectory u(t) stabilizes as t — » oo to the set of equilibria 

K := {u £ H 2 , Au - Uf(u ) = Tig}. (3.8) 

Furthermore, if the setlZ is discrete, every trajectory u(t) converges to a single equilibrium 
uq £ TZ and the rate of convergence is exponential if that equilibrium is hyperbolic. 

Remark 3.8. Note that, for generic g £ L 2 , the set 1Z will contain only hyperbolic equi- 
libria (see PP ) . In that CcLSG, clS it is not difficult to prove (again verifying the conditions of 
the abstract theorem on regular attractors stated in p]), the attractor A can be presented 
as a finite union of finite-dimensional submanifolds of H (the unstable manifolds of all 
equilibria) and that the rate of attraction of any bounded subset B to the global attractor 
A is exponential. 

4. The convective Brinkman-Forchheimer equations 

In this section, we extend the results of the previous section to the case of the following 
Brinkman-Forchheimer equation with the Navier-Stokes type inertial term: 

d t u+(u,V)u + Vp = Au- f(u)+g, divw = 0. (4.1) 

Note that the case / = corresponds to the classical Navier-Stokes problem and the 
general case / ^ can be also considered as the so-called tamed Navier-Stokes equation, 
see P2|. 

As before, the nonlinearity / is assumed to satisfy conditions (II. 3p but with the ad- 
ditional lower bound r > 3 which is necessary for the uniqueness. Note that no upper 
bounds for the growth exponent is posed. 

As before, we define a weak solution u as a function of the class (13. ip satisfying (14.11) in 
the sense of distributions, see Definition 13. II In addition the assumption r > 3 guarantees 
that 

(u, V)u £ L 4/3 C L q , q:={r + 1)* < 4/3 (4.2) 

and, therefore, in contrast to the case of the classical Navier-Stokes equations, the mul- 
tiplication of (14. ip by u with integration over Q is justified for any weak energy solution 
of that equation. Thus, we have verified that any weak energy solution of (14.11) satisfies 
the energy estimate (12.11) . The existence of an energy solution can be then obtained in a 
standard way via the Galerkin approximation method. 

The next Lemma gives the uniqueness of the energy solution for the case r > 3. 
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Lemma 4.1. Let the nonlinearity f satisfy (jl.3p with r > 3 and g G L? . Then, for every 
uq G H , the problem (14. ip possesses a unique weak solution u and this solution satisfies 
the energy estimate (12.11) . 

Proof. Indeed, let U\ and u 2 be two solutions and let v — U\ — u 2 . Then, this function 
solves 

d t v + {v,V)u 1 + {u 2 ,V)v + Vq = Av-[f{u l )-f{u 2 )}, divi; = 0. (4.3) 
Multiplying this equation by v, integrating by parts and using that / satisfies (11.31) . we 
will have 

j t \\v\\ 2 L2 + 2\\Vv\\l 2 + adu.r 1 + \u 2 r\ \v\ 2 ) < C\\v\\% + 2\((v, V)m,v)\ 

for some positive a depending on k from ( II. 3ft . Here we have implicitly used that the first 
condition of ( 11 .3p implies that 

if{ui) - /(t*a),ui ~u 2 ) > -C\\u! -u 2 \\ 2 L 2 + a(|wi| r_1 + |u 2 r~\ |«i ~u 2 \ 2 ), 

see p/4] and [5] for the details. 

The last term in the above differential inequality can be estimated integrating by parts 
once more and using that r — 1 > 2: 

2\{(v, V)u h v)\ < 2{\ Ul \ ■ \v\, \Vv\) < \\Vv\\ 2 L2 + C(\ Ul \ 2 , \v\ 2 ) < 

< \\Vv\\ 2 L2 + a(|u 1 | p - 1 + |u 2 r -1 , M 2 ) + C\\v\\ 2 L2 . (4.4) 

Thus, we have 

||HI| 2 + I|V^|| 2 , 2 <C'|H| 2 L2 (4.5) 

and the uniqueness is proved. □ 

Remark 4.2. As we see from the proof, the uniqueness holds for the case r = 3 if 
the coefficient k in ( II .3p is large enough. However, we do not know whether or not the 
uniqueness holds for any cubic nonlinearity (without this assumption). 

The next theorem is analogous to Theorem 13.31 and gives the regularity of solutions for 
problem (I4.ip . 

Theorem 4.3. Let the function f satisfy ( 11. 3p with r > 3 and let g G L 2 . Then, for any 
Uq G H, the associated solution u(t) of (14. ip is more regular for t > (u(t) G H 2 ) and 
estimate (12. 16f) holds. In addition, if uq G H 2 then estimate (12. 151) also holds. 

Proof. The proof of this theorem is also analogous to the proof of Theorem 13.31 Indeed, 
differentiating equation ( 14.11) with respect to t and arguing as in the proof of the previous 
lemma, we conclude that the function v = dtu satisfies the differential inequality (14.51) . 
On the other hand, using (14. 2 p for the control of the inertial term and arguing as in 
Corollary 12. 3^ we derive estimate (12.31) and based on that estimate and inequality (14.51) 
for v = dtu, one derives the controls (12. lip and (12. 7p for the time derivative v = dtu (all 
these estimates can be justified via the Galerkin approximations). 

Finally, having the control of the L 2 -norm of d t u, one can treat problem (14.11) as an 
elliptic boundary value problem of the form (I5.42p and apply Corollary 15.41 which gives 
the desired estimate for the if 2 -norm and finishes the proof of the theorem. □ 
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Remark 4.4. Note that the nonlinear localization technique used in the proof of Corollary 
15.41 is not necessary if r < 5 where we may use the standard maximal regularity for the 
linear Stokes equation or in the case of periodic boundary condition. However, we do not 
know how to avoid these technicalities in a general case. 

Finally, let us note that the analogue of Theorem 13.51 holds for the Navier-Stokes case 
as well. 



Theorem 4.5. Let the assumptions of Theorem \4-^\ hold. Then the solution semigroup 



S(t) : H —¥ H possesses a global attractor A which is a bounded subset of H 2 and possesses 
the standard description (13. 6p . 

The proof of this theorem repeats word by word the proof of Theorem 13.51 and so is 
omitted. 

Remark 4.6. To conclude, we note that all assertions formulated in Remark 13.61 remain 
true for the convective case as well. 

5. Appendix: Maximal regularity for semi-linear Stokes problem 

The appendix is devoted to the stationary problem associated with the problem (11. ip . 
that is the following semi-linear Stokes problem: 

-Aw + f(w) + Vp = g, divw — 0, x G f2, ^ ^ 

w = 0, x G dQ, J n p(x) dx = 0. 

Here w = (wi,w 2 ,w 3 ), g G L 2 (Q) is a given function and the nonlinearity / satisfies 
assumptions (II .3p with arbitrary r > 1 and with K — 0. Thus, we have assumed that 
the nonlinearity / is monotone f'{u) > and, therefore, the energy solution of (15.11) is 
unique. 

Our aim here is to prove the L 2 -maximal regularity estimate for problem (15.11) (which 
is the non-linear version of the classical L 2 -regularity theorem for the Stokes operator). 
Before stating the main result, we first remind the straightforward L 9 -regularity result 
where q = (r + 1)* = 1 + K 

Lemma 5.1. Let the above assumptions on f hold and let g G L q (Q). Then, problem 
(15. ip has a unique solution (w,p) G T q where 

jr ;= {( w , p ), w e W 2 > q ({l) n L r+ \Q), p G W 1,q (Q)} 

and the following estimate holds: 

\\w\\ W 2, q + ||to||ir+i + IMIl3/2+ £ + ||w|| 2 ^ 1 /(r+1) < C(l + \\g\\Li), (5.2) 

for some positive C independent of g and sufficiently small e = e(q) > 0. 

Proof. We give below only the derivation of the estimate in the space T q (the existence 
and uniqueness of the solution can be obtained in a standard way, e.g., using the Galerkin 
approximation method). Indeed, multiplying equation (15. ip by w, integrating by parts 
and using (II. 3p . we arrive at 

IMIfli + \H\ r itl < C(l + \\g\\l) , q:=(r + iy = l + i (5.3) 
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where C is independent of g and w. Together with conditions on /, this gives, in particular, 
that 

\\f(w)\\ Lq <C(l+\\g\\ Lq ), (5.4) 

for some (new) constant C. 

Rewriting now the problem (15 .ip as a linear Stokes problem 

- Aw + Vp = h w := g - f{w) (5.5) 

and applying the maximal L 9 -regularity estimate for this linear Stokes problem, we have 

||w||vk 2 .9 + IMIvy 1 ^ < C(l + IMUO- (5-6) 

In particular, due to Sobolev embedding theorem W 1,q (£l) C L s (£l) with s := > 3/2 
(since r > 1 and q < 1) 

IHIl3/2 +£ <C(1 + |M| l ,), (5-7) 

where e = e(r) > depends only on the exponent r. □ 

We are now ready to state the main result of this section. 

Theorem 5.2. Let w be an energy solution of problem (15. 1 j) . g € L 2 (Q) and the assump- 
tions (11.31) on f hold. Then, w G H 2 (Q) and the following estimate is valid: 

IMMfi) + \\p\\m(n) < Q(\\(w,p)y q )(l + ||^||^ K ) (5.8) 
for some monotone function Q and positive k = k{t). 

Proof. As before, we restrict ourselves to the formal derivation of the regularity estimate 
(15. 8p . The existence of a solution can be verified in a standard way using, e.g., the 
Leray-Schauder fixed point theorem 

We will use the so-called nonlinear localization method and split the derivation of that 
estimate in several steps. 
Step 1: Interior regularity. 

At this stage, we obtain the interior if 2 -regularity estimate for the solution w. To 
this end, we multiply equation (15. ip by d Xi (4>d Xi w) where is a proper nonnegative 
cut-off function which vanishes near the boundary and equals one identically inside of the 
domain. To be more precise, we assume that if G C^°(IR 3 ) is such that < <f(x) < 1, 
(p(x) = if x G M 3 \f2 and f{x) = 1 if x G Q u where 

Vl u := {x G fi, dist(x, d£l) > v} 

and v > is a sufficiently small number. Moreover, without loss of generality, we may 
assume that 

|V0(x)| < C^Or) 1 - 5 (5.9) 

for some 5 > which can be chosen arbitrarily small and the constant C U) s depending 
only on 5, v and the shape of Q. 
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Thus, we may estimate the term with the Laplacian as follows: 
(Aw,^2d Xi (tpd Xi w)) = ^2(Vd Xi w,V((pd Xi w)) > K^2(ipVd x .w,Vd Xi w)- 

i i i 

- {Vd Xi w, V<p ■ d Xi w) > ^K{<p, \Dlw\ 2 ) - (|V^| V 1 , |Vw| 2 ) > 

> ^y^wWlz -c\\wf m . (5.10) 

This, together with the energy estimate (15.31) for the subordinated terms, gives the fol- 
lowing estimate: 

U ll2 w\\ 2 H* + (<t>f'(w)Vw,Vw) < C(l + \\g\\ 2 L2 ) + \(p,d Xi (V(f>-d Xi w)\. (5.11) 

Using again the energy estimate to control the subordinated terms and (15. 9p to control 
\Vif\, the last term can be estimated as follows: 

|(p, <9, S (V0 • d Xi w)\ < C(l + \\g\\l 2 ) + eU^wWl* + ^||0 1/2 ^||i 2 , (5.12) 

where e > is arbitrary. In order to estimate the last term in the right-hand side of this 
estimate, we use the Holder inequality with exponents ^ and (a = 1/2 — S) in the 
following way: 

/ v 2a M 2 dx = [ ( V \u\) 2a \u\^- a Ux < Cy U f?4uf%%. (5.13) 
Jq Jq l 

Since ^? = |(1 + ^) and 2a = 1 - 25 < 1, fixing 5 > small enough that ^ < §£ 
and using ( 15. 71) for estimating the L 3 / 2+£ -norm of p, we arrive at 

H 1/2 - 5 P \\h < c\\M\inp\\t^ < 0(IIKp)lk)(i + \\M»)- (b.h) 

Thus, due to ( 15TT2|) . 

U 1/2 M\hi + WMVw, Vw) < C(l + |b||| 2 ) + 3(11(^)11^X1 + ||#|Ua). (5.15) 

In order to estimate the term in the right-hand side, we take the divergence from both 
sides of ( 15. ip and write out 

Ap = — div/(w) + divg. (5.16) 
Multiplying this equation by 0, we have 

A(4>p) = 2V0 • Vp + A(f>p - (f)divf(w) -(pdivg := h. (5.17) 
Furthermore, due the growth assumptions (II. 3p on / and the energy estimate (15. 3p . 

\\<Pdivf(w)\\l q < [ (j) q \f{w)Vw\ q dx< [ \<f)f(w)Vw-Vw\ q/2 -\<i)f(w)\ q/2 dx< 
Jq Jq 

< C{cf ) r{w)Vw^w) q l 2 \\f{w)f\_ = C(<j>f(w)Vw,Vwy/ 2 \\f(w)\\ q/2 +1 < 

<Q(\\(w,p)\\^)(l + (cl>f(w)Vw,Vw)y/ 2 . (5.18) 
Using this estimate together with the energy estimate for the pressure, we arrive at 

WHu+h-i < /3(0/'HV«;,V«;) + g^||(«;,p)|U)(l + |b|U 2 ), (5.19) 
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where the positive constant may be chosen arbitrarily small (the term H^ 1 appears due 
to the term (pdiv g, other terms belong to L q ). Finally, due to the maximal regularity for 
the Laplacian together with the Sobolev embedding W 2 ' q C L s with s = t^tt > 3 and 
H 1 C L 6 , we have 

\\M\w^+m + ||#||L3+e < P(<j>f'(w)Vw,Vw) + Q p (\\(w,P\\r<)(l + IMM, (5-20) 

for some small positive e = e(r) and 0. Inserting this estimate into (15.151) and fixing j3 
small enough, we conclude that 

U 1/2 w\\h + W(w)Vw,Vw) < Q(\\(w,p)\\s q )(l + 1Mb)- (5.21) 

Using the embedding H 2 C C and the fact that <fr equals one identically inside of the 
domain, we deduce the desired interior regularity 

IM|fl*(n„) + \\p\\m(n v ) < Qu(\\(w,p)\\^ q )(l + ||^||^a), (5.22) 

where v > is arbitrary. Thus, the interiror if 2 -regularity estimate is proved. 
Step 2: Boundary regularity: tangent directions. 

We now obtain the if 2 -regularity in tangent directions near the boundary. The standard 
approach here is to make the change of variables which straighten the boundary locally 
in a small neighborhood of a boundary point xq and then obtain the global estimate 
using the proper partition of unity. However, in order to avoid the complicated notations, 
we will use the alternative equivalent approach working directly with the derivatives in 
tangential directions. Namely, let T\ = r x (x) and r 2 = r 2 (x) be two smooth vector fields 
in Q which are linear independent in a small neighborhood Q\Q U and such that, for any 
Xq G dQ u , v G [0, u ] the vectors Ti(xq) and r 2 (x ) generate the tangent plane to dQ u . 

Being pedantic, such vector fields usually do not exist globally, but only locally (in a 
neighborhood of a fixed point xq G dQ). However, the plane-field spanned by the pair 
of vectors (n, r 2 ) is well-defined globally and the tangent gradient V T is also well-defined 
globally. Nevertheless, in slight abuse of rigoricity and in order to avoid the completely 
standard technicalities, we assume that the both vector fields T\ and r 2 are globally defined 
in Q. 

Let z(x) := d T w(x) := Y^H=x r% { x )^ i w (where r = T\ or r = r 2 ). Then, this function 
solves 

-Az + f'{w)z + V(d T p) = h z , 

h z := d T g + d Xi (Ti(x)p + Ki{x)Vw) + LVw + M(x)p, (5.23) 

div z = C{x)Vw 1 z\ m = 

for some smooth (matrix) functions T, K, L, M, C independent of p and w. Multiplying 
this equation by z and using the energy estimates (15. 6p and (15.71) (in order to estimate 
the subordinated terms) together with the facts that z\ m = and div z = C(x)Vw, we 
deduce after the simple estimates that 

\\z\\ 2 Hl + (f'(w)z,z) <Q(||(^p)||^(l + |M|| 2 ) + C'||V W |U3||^|| X 3/ 2 + C'|b|| 2 , 2 . (5.24) 

Let us first estimate the L 3 -norm of Vw. To this end, we need to use the interpolation 
in the spaces with different regularity in tangent (t± and r 2 ) and normal (n) directions 
(here we need also the boundary dQ to be smooth enough). Indeed, since we control the 
W 2,q -norm of w (due to estimate (15.61) ). from the Sobolev's trace theorem, we have the 
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control of the traces of Vw on the surfaces dVl u in the H^ 1 / r ' 9 -norm (remind that for small 
v, dVLy are uniformly smooth if the initial dVt = <9f2 is smooth and the Sobolev trace 
theorem as well as the interpolation theorems below work): 

<Q(\\(w, P )h q ). (5.25) 

On the other hand, using the estimate (15.241) for Vd T w, and the embedding H 1 C L s for 
every s < oo for 2D domains Q u , we conclude that 

II Vw\\L2(v£[o,v ],L°(dn v )) < Cspllifi, (5.26) 

where the constant C s depends only on s. 
Using now the standard interpolation 

(L 1+1 / r (dn u ),L s (dn u )) 2/3 = L 3 (dn u ) 

if | = | i+i/r + Is' *- e -' w ith s r = 2(1 + r), we see that 

||Vu>|| L 3 ( n\fi„ ) = IIVwIIls^jo^^s^)) < 

< C r ||VwI|^ (ver0j ^ ] ^ 1+1/ ^ ( ^ )) ||Vw||^ e[0j ^ ]f ^ (oni , )) < Q{\\{w,p)\\r t )\\4 2 Hl (5-27) 
Using also the interior regularity (I5.22p . we conclude that 

\\Vw\\ LSm < Q(\\(w,p)\\^)(l + \\z\\f l{n) ). (5.28) 

Let us now estimate the L 3//2 -norm of d T p. To this end, we rewrite equation (15.231) in the 
form 

Az - V{d T p) = h, z\ 9n = 0, divz = C(x)Vw (5.29) 

with 

h = f{w)z + d T g + d Xi {Ti{x)p + Ki(x)Vw) + L(x)Vw + M(x)p (5.30) 
and note that, due to the energy estimates (15.61) . (15.71) and similar to (I5.18p . 

k-i+w < Q(HKp)|k)(l + (f'(w)z,z) 1 / 2 ) + C\\g\\v. (5.31) 



Applying the H 1 — > H 1 regularity and L q — > W 2,q regularity estimates for the linear 
non- homogeneous Stokes equation, similar to (I5.20p . we conclude that 

\\d T p\\ L 2 +w x, q + \\d T p\\ L3/2+ , < Q(\\(w,p)\\ Tq )(l + (f(w)z, z) 1 ' 2 ) + C\\g\\ L 2, (5.32) 

where e = e(r) > is small enough. 

Finally, we need to estimate the L 2 -norm of p. To this end, using the fact that we have 
the control of the W^^-norm of the pressure p and Sobolev trace theorems, similar to 
(15.25!) . we have 

\\p\\L™{ve[o, Vo ],Li(da„)) < Q(\\(w,p)\\ Tq ). (5.33) 
On the other hand, due to the embedding W 1 ' 3 ^ 2 C L 6 for the 2D domains Q u , 

IblUs^elo^ol.^n,)) < C\\d T p\\ L3 / 2 (5.34) 
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and, consequently, due to interpolation with exponent (L 9 / 4 = (L 1 , L 6 ) 2 /s = (L°°, L 3 / 2 ) 2 /3) 
together with already proved interior regularity (and the fact that 9/4 > 2), we deduce 
the estimate: 

L 2 (Q) < |blU 2 (fi\O, ) + lbl|L 2 (n„ ) < 

<Q(IIKP)lk)(l+ll 



< Q(\\(w,p)\\ Tq )(l+\\g\\ L 2 + |b|li / i (l/eMiL i ( an I/)) lbll^ 3 / a(l , 6[ o^ ]iL 6 (a n„))) 



< 



<Q(\\(w,p)\\ Tq )(l + \\g\\ L , + \\d T p\\%%) 



< 



< Q(\\(w,p)y q )(l + \\g\\ L 2 + (f'(w)z,z) 1 / 3 ). (5.35) 

Inserting now estimates (15. 28ft . ( 15.321) and (I5.35P into the right-hand side of ( (15. 24ft . we 
finally arrive at 

\\4m + HVHIia + (f'(w)z,z) + |b||| 2 < Q(\\(w,p)\\ Tq )(l + |b||| 2 ) (5.36) 

and the .^-regularity of w in tangent directions is verified. Now, we note that 

\\f(w)z\\ H -i < C\\f'(w)z\\ L6/5 < C(\f'(w)\z,z))^\\f'(w)\\f /4 

and, therefore, applying the if -1 -regularity theorem to the linear Stokes problem (15.291) 
and using (15.361) . we arrive at 

\\d T p\\v < Q(IIKp)lk)(i + \\g\\v)(i + \\f{w)\\%). (5.37) 

Step 3: Regularity in normal direction and the final estimate. Let us now 

multiply equation (15. ip by Aw. Then, after integration by parts, we get 

IMIh* + (l/'HlVw, Vm) < C(l + \\g\\ 2 ) + C\(Vp, Aw)\. (5.38) 
In addition, the second term in the right-hand side gives 

(f'(w)Vw,Vw) > K(\w\ r -\ \Vw\ 2 ) > /«i||V(|«;| (r+1)/2 )||i 2 > 

>*HI3ki> >K4f(w)& +1) -c 

and, consequently, using the energy estimate and the interpolation 

||/HI| i2 <c7||/H||^||/H||^ 4 , 

we will have 

II* + H/MllS r+1)/(3r) < C(l + \\g\\ 2 ) + C\(Vp, Aw)\. (5.39) 

Thus, we only need to estimate the last term in the right-hand side of this inequality. To 
this end, we split the tangential and normal derivatives in that term and use (15.361) and 
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(I5.37P for estimating the tangential derivatives: 
\{Vp,Aw)\ < \{d T p, (Aw) T )\ + \{d nP , {Aw) n )\ < ^\\w\\ 2 H2 + 

+ C\\d T p\\ 2 L2 + \(d nP , d 2 n w n )\ + C(\V P \, \Vd T w\ + \Vw\) < ~\\w\\ 2 H2 + 

+ Q(\\(w,p)\\^ q )[l + (1 + \\gf L2 )\\f(w)\\ L5/4 + ||Vp|U 2 (l + \\g\\ L >)}+ 

+ C\\V P \\ L 4d 2 n w n \\ L 2. (5.40) 

To estimate the last term in the right-hand side of (15.401) . we use that w is divergent free 
and, therefore, 

d n w n + d Tl w Tl + d T . 2 w T2 = C(x)w. 

Differentiating that equation in the direction of the normal vector field, we obtain the 
estimate 

\\dlw n \\ L 2 < C\\Vd T w\\ L 2 + C\\Vw\\ L 2. 

Inserting that estimate to the right-hand side of (15.401) and (I5.39j) and using (I5.36P together 
with the obvious estimate 

\\Vp\\ L 2<C\\f(w)\\ L 2+C\\g\\ L 2 

(which follows from the L 2 -maximal regularity for the linear Stokes equation) and the 
interpolation 

ll/HILvs<c||/H||^ 5 ||/H||^, 

we arrive at 

\H\h + \\fM\\T m3r) < 

< Q(ilKp)lk)[i + (i + \\g\\h)\\fW\\%* + WfWWvO- + bWv)]- (5-4i) 

Finally, thanks to Young inequality, we derive from (I5.4ip that 

\\w\\ H2 <Q(\\(w,p)\\ Tq )(l + \\g\\ 2 L r), 
where k = /t(r) > is a positive number. Thus, the theorem is proved. □ 

Remark 5.3. Clearly, when the nonlinear terms has the growth 

\f(u)\ < C(l + |m| 3 ),Vm G M 3 , 

the maximal regularity estimate (I5.8P follows directly from the energy estimate (15.31) and 
the regularity estimate for the linear Stokes problem. Note also that, in the case of periodic 
boundary conditions the simple multiplication of the initial problem by Aw gives better 
estimate 

\\W\\ H 2 < C(l + ||#|| L 2). 

However, for the case of Dirichlet boundary conditions, the additional (uncontrollable) 
boundary terms appear under the integration by parts, and we unable to obtain the ir- 
regularity estimate which is linear with respect to the L 2 -norm of g. However, as we will 
see below, the sub-quadratic growth rate of that estimate with respect to g is enough to 
be able to apply it for the Navier-Stokes-type problem. 
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To be more precise, we want to apply the above result to the following analogue of 
problem (15. ip perturbed by the Navier-Stokes inertial term: 

-Aw + (w, V)w + f(w) + Vp = g, div w = 0, x G f2, ^ . 

w = 0, x G dQ, J n p(x) dx = 0. 

Indeed, since the inertial term vanishes after the multiplication the equation by w and 
integrating over x, arguing as in Lemma 15.11 we have: 

HI^ + IMI^^i + NIU 

In addition, 

\\(w,V)w\\ L3/ 2 < C\\w\\ L6 \\Vw\\ L 2 < Q(\\g\\v), (5.43) 

for some monotone function Q. Thus, the L 9 -norm of the inertial term is under the 
control if q < 3/2 (= r > 2) and, applying the /^-regularity estimate for the linear Stokes 
problem, we also have the control of the W 1,q -norm of p. So, we have proved that 

\\(w,p)U<Q(\\9\\w) ( 5 - 44 ) 
if r > 2. The next Corollary gives the if 2 -regularity estimate for the problemf l5.42p . 

Corollary 5.4. Let f satisfy (jl.3p with r > 2 and g G L 2 . Then, any energy solution 
(w,p) G of problem (I5.42p belongs to H 2 x H 1 and the following estimate holds: 

IMIaa + Ibll^ 1 < <2(IMU 2 ) ( 5 - 45 ) 
for some monotone function Q. 

Proof. Let us first formally deduce a priori estimate (15.451) . To this end, we interpret the 
inertial term as an external force and apply estimate (15.81) . Then, using (15.441) . we have 

\\w\\ H 2 <Q(\\g\\ L ,)(l + \\( W y)w\\ 2 L ^) (5.46) 

for some monotone Q and some positive k = n(r). Thus, we only need to estimate the 
L 2 -norm of the inertial term. To this end, we use (I5.43P together with the interpolation 
inequalities and the fact that H 2 C W l,& : 

\\(w,V)w\\ L . < \\(w,V)w\\f /2 \\(w,V)w\\)! b 3 < 

< QMrfM&WfS < Qi(h\\L*)\H\)iUM\¥i < Q2{\\g\\»)M l £- ( 5 - 47 ) 

Inserting this estimate in the right-hand side of (15.461) . we deduce the desired a priori 
estimate (15.451) . 

The existence of a solution (w,p) G H 2 x H l can be obtained in a standard way based 
on that estimate and approximating, for instance, the growing non-linearity f(w) by a 
sequence f n (w) of globally bounded ones (see, e.g., [9] for the details). However, since the 
solution of (15.421) may be not unique, we still need to verify that any energy solution of 
that equation satisfies estimate (I5.45p . 

Indeed, let w be an energy solution of (I5.42p . Let us consider the following modified 
equation (15.421) : 

-Av + (v,V)v + f(v) + Rv + Vp = g w , divv = 0, x G fi, 4 ^ 
v = 0, x G dQ, J n p(x) dx = 0. 



ATTRACTOR FOR THE BRINKMAN-FORCHHEIMER EQUATIONS 



19 



with g w := g + Rw. We claim that the solution v — w of that equation is unique in the 
class of energy solutions if R 3> 1 is large enough (this fact can be easily verified using 
the standard energy estimates). On the other hand, arguing as before, we can construct 
a regular solution (v,p) G H 2 x H 1 of that equation satisfying ( 15. 45ft . The uniqueness 
guarantees then that the initial solution u is also regular and satisfies this estimate. Thus, 
the corollary is proved. □ 

Remark 5.5. Of course, the condition r > 2 is not necessary for the validity of Corollary 
15.41 Indeed, for r < 2, the L 2 (and even L 3 ) norm of nonlinearity f(w) is under the 
control due to the initial energy estimate and the desired if 2 -regularity can be obtained 
treating the nonlinearity / as the external force (exactly as in the case of the classical 
stationary Navier-Stokes problem, see e.g. [I0])- However, we are mainly interested by 
the case of large q where the above presented nonlinear localization technique becomes 
unavoidable. 
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